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Abstract 

Coherent electron transport in double-barrier heterostructures with parallel 
electric and magnetic fields is analyzed theoretically and with the aid of a 
quantum simulator accounting for 3-dimensional transport effects. The onset- 
voltage shift induced by the magnetic field in resonant tunneling diodes, which 
was previously attributed to the cyclotron frequency w c inside the well is 
found to arise from an upward shift of the non-zero ground (lowest) Landau 
state energy in the entire quantum region where coherent transport takes 
place. The spatial dependence of the cyclotron frequency is accounted for and 
verified to have a negligible impact on resonant tunneling for the device and 
magnetic field strength considered. A correction term for the onset-voltage 
shift arising from the magnetic field dependence of the chemical potential is 
also derived. The Landau ground state with its nonvanishing finite harmonic 
oscillator energy hw c /2 is verified however to be the principal contributor to 

the onset voltage shift at low temperatures. 
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I. INTRODUCTION 



One of the most fundamental results of quantum theory is the fact that the harmonic 
oscillator has a ground state with a non-zero energy hw c /2. This energy is usually ignored 
in most field calculations because it is merely a constant and it is therefore treated as a 
reference point in energy. 

The quantum state of an electron in a constant and uniform magnetic field is described 
by a harmonic oscillator with cyclotron frequency w c . Indeed from a classical mechanics 
point of view, this results from the fact that the electrons move through the device according 
to helicoidal trajectories, and that the projected motion along an arbitrary transverse axis 
is that of a simple harmonic oscillator. Varying the magnetic field changes the cyclotron 
frequency and is found in turns to induce a shift of the onset voltage in the IV character- 
istic of double barrier heterostructures. We demonstrate in this paper that this collective 
phenomenon in the electron current occurs due to the varying energy value hw c /2 of the 
ground-state of the harmonic-oscillator representing the electrons (see Fig. [3]). 



magnetic field is represented by a harmonic oscillator. The finite ground state energy of an 
electromagnetic field between parallel metal plates depends on the separation of the plates, 
which is similar to the dependence upon the varying magnetic field in the system studied 
here. The Casimir effect which was predicted in 1948 was recently verified experimentally 
by Lamoreaux in 1996@. 

Other known examples are the Lamb shift!, and the unfreezing of liquid helium at ex- 
tremely low temperatures. Nevertheless phenomena which have a direct relation with the 
nonvanishing, finite ground-state (zero-point) energy of the corresponding harmonic oscilla- 
tor are very few. In this paper one more occurrence is demonstrated for magnetotransport 
in semiconductor heterostructures. 

The magnetotransport phenomenon in multiple barrier heterostructures with a magnetic 
field applied parallel to the electric field was first investigated and reported in the pioneer- 
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same origin in nature. Each mode of the electro- 
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ing work of Mendez, Esaki and Wangu. When both the magnetic field and electric field 
are perpendicular to the heterostructure interface!, the electron motion becomes quasi-one 
dimensional^, as the electron states in the plane perpendicular to the layered direction are 
quantized in discrete Landau levels. This dimensional reduction of the electron motion has 
a profound effect on the transmission probability and current- volt age characteristic when 
compared with the case of zero magnetic field. 

The experimental work on GaAs/GaAlAs double barrier structures incorporating both 
a parallel electric field and a magnetic field showed some remarkable features, including 
a shift of onset voltage and enhanced peak-to-valley-ratic&i. However, these observations 
have not been fully accounted for from a theoretical standpoint, and such a theoretical 
analysis is needed to understand the underlying physical processes and inspire possible 
device applications. 

Specifically in Ref. 4 the shift of onset voltage with increasing magnetic fields was 
predicted to be given by the energy hw c /2 associated with the cyclotron frequency inside the 
well region. In this paper we demonstrate instead that the onset voltage shift originates from 
an upward shift of the non-zero ground (lowest) Landau state energy in the perpendicular 
plane throughout the entire region of the quantum structure where coherent transport takes 
place, including the left contact, rather than only inside the well. We also derive and 
numerically calculate a correction term for the onset voltage shift which arises from the 
magnetic field dependence of the chemical potential. Finally we have verified using our 
3D transport model that the spatial dependence of the cyclotron frequency has a negligible 
impact on resonant tunneling in the test device for the magnetic field strength considered. 

From a fundamental physical point of view, the fact that the nonvanishing ground state 
energy 1/2 • hw c reveals itself in a shift of the onset voltage is a very noticeable phenomenon. 
We account here for this abnormal feature with the aid of a quantum simulator based 
on a 3-dimensional transport picture which incorporates the Landau states in the 
perpendicular plane. 

Undoped Alo.3Gao.7As/GaAs/Alo.3Gao.7As double barrier structures with heavily doped 



n + -GaAs left and right contact layers are considered. The donor doping density is 1 x 
10 18 /cm 3 . The conduction-band profile without an applied voltage is shown in the inset 
of Fig. [I]. The effective mass of the electrons for GaAs and Alo.3Gao.7As are given by 
m* = 0.067m, 0.092m, respectively, where m is the electron rest-mass. 

Two structures with a barrier and a well width of 34A (structure A) and 50A (structure 
B) respectively are examined. We assume that the donors in the contacts are completely 
ionized (n = No) due to the heavy doping density. A tight-binding band in the superlattice 
direction and a parabolic band in the perpendicular plane are used. Spin splitting in a 
magnetic field is ignored for simplicity. Finally we focus on the low temperature case, 
T = A.2K. 

II. THEORY 

The electron Hamiltonian in the heterostructure in the most general case includes inter- 
action terms which describe disorder scattering by phonons, impurities, interface roughness, 
etc. Since the electron mean free path at very low temperatures is relatively large compared 
to the heterostructure size under investigation (w c t ^> 1), we can safely neglect scatter- 
ing. Therefore, only the coherent, ballistic electron motion will be considered here and the 
Schrodinger equation to be solved is 

H e \*)=ih±\*). (1) 

In the plane perpendicular to the superlattice direction, the electrons move classically in 
a circular path. In a quantum treatment the electron energy is quantized in Landau states!, 
which are given by the quantized harmonic oscillator states, 

E p = (L + l/2)hw c (n), (2) 

where w c (n) = eB/m*(n) is the cyclotron frequency in the magnetic field B at the lattice 
site n in the superlattice direction. In the following discussion we will see that the lowest 
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Landau state L = is not just a reference point but gives an important contribution to the 
current flow in the double barrier heterostructure. 

We expand the electron wave function in a basis consisting of generalized Wannier func- 
tions in the superlattice direction (x-direction)il and of the Landau states in the (y, z) plane 
which result when the Landau gauge is chosen for the vector potential of the magnetic field, 

l*> =Y, j f( k yi L ' n )\ k yi L i n ) dk yi ( 3 ) 

n.L 

where f(k y , L, n) is the envelope function of the electron states. The matrix element for the 
heterostructure Hamiltonian between generalized Wannier-Landau states is§ 

(k y , L, n\H e \k' y , L', n') = H e , nn ,(L)5 L , L i5(k y - k' y ), (4) 

where 

H e)nn i(L) = H nn > + (L + l/2)%w c (n)6 nn >, (5) 



and 



H * 

2a 2 J m* {n)m* (n + 1) 



h 2 



2 



+ E con (n) - eV am (n) 



o \ 1 ^cony ^ r app\ 



<W- (6) 



2a 2 y ^m* {n)m* (n — 1) 

Here E con {n) is the conduction band edge and V app (n) is the electrostatic potential at the 
lattice site n. 

The transversal mass in the perpendicular part of the Hamiltonian matrix also varies 
along the longitudinal direction. We introduce a new Hamiltonian matrix which absorbs 
this n-dependent transversal part as 

#e,rm'(£) = H nn ,{L) + (L + l/2)7iW c (0)<W, (7) 

where 

H nn >(L) = H nn , + (L + l/2)hw c (0)[m*(0)/m*(n) - 1]<W- (8) 



Here H nn i depends on the perpendicular states L, incorporating 3-dimensional band- 
structure effects in the longitudinal motion of the electrons. 

Thus the Schrodinger equation reduces to the following 3-dimensional envelope equation 
for f(k y , L,n): 

m df(k y L,n) = (L + l/2)hWcof(k ^ L , n ) + £ H nm f(k y , L, m), (9) 
at m 

where w c0 = eB/m*(0) is the cyclotron frequency at the left contact. 

The incident electrons in the left contact occupy the sub-bands E Lo = E 0x + (L + 
l/2)hw c0 , where E 0x = h 2 kl/(2m*(0)). Since in the absence of scattering an incident ballistic 
electron conserves its Landau level L and its total energy E Lo , as it goes through the 
structure, the envelope function solution of Eq.(9) in Landau gauge takes the form 

f(k y , L, n) = e-^o* f (n)5(k y - k y0 )6 LLo . (10) 

When Eq.(10) is substituted into the Schrodinger equation (9), we arrive at the time inde- 
pendent form 

E 0x fo(n) = H n mfo(m) (11) 

m 

Note that the longitudinal electron wave function fo(n) depends on the perpendicular Lan- 
dau sub-band level L , since the second term in the right hand side of Eq.(8) is a function of 
the Landau level. This is a 3-dimensional transport effect which was exactly implemented 
in our numerical analysis. 

The electron currents for an incident and transmitted state (E 0x ,L ) are respectively, 

j I0 (E 0x ,L ) = ev L (E 0x ) 

3T0 (E 0x ,L ) = e\f {n R )\ 2 v R (E 0x ), (12) 

where vl,vr are the electron velocities at the left and right contacts, and tlr is the lattice- 
site index of the right contact. The transmission probability from left to right contact is 
then defined as 
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Tl^r(E 0x , Lq) = jro(Eox7 L )/j I0 (E 0x , L ). (13) 

It is noted that the longitudinal transmission probability depends on the perpendicular state 
Lq. Fig. |I] shows the transmission probability for different applied biases. 

In order to calculate the total diode current, the chemical potential corresponding to 
a doping density n in the contacts in the presence of a magnetic field should be deter- 
mined. The chemical potential has an oscillatory dependence on the magnetic field which is 
determined from the carrier density through the following integration relation, 

n = N d = {V2^*/irh)D £ / (1/ E 0x ) f FD (E 0x , L , ft(B))dE 0x , (14) 

L =0 J ° 

where D = eB/h is a degeneracy factor in the plane perpendicular to the superlattice 
direction and where the Fermi-Dirac function f FF > is given by 

f FD (E 0x ,L ,fi) = {exp[(E 0x + (L + l/2)hw c0 - n)/kT\ + l}" 1 . (15) 

The chemical potential /i is a complex function of the electron density n, magnetic field and 
temperature. Its numerical calculation shows that the chemical potential [i oscillates by a 
large quantity when the magnetic field B is strong; for example, there is a difference of 4.4 
meV between its value at B = 13 and at B = 17 Telsa. It is important to account for this 
behavior to obtain an accurate value of the current. 

Under strong magnetic fields, the density of electronic states (see Eq.(14)) drastically 
departs from the magnetic-field free case and in addition an upward shift hw c (n)/2 of elec- 
tron's lateral energy with respect to the conduction band minimum occurs in the entire 
quantum region. As we will see below, the current integration over electronic states in 
the local emitter will strongly depend on the Tiw c q/2 shift via the Fermi-Dirac distribution 
function, Eq.(15). This will lead in turn to the onset voltage shift in the current-voltage 
curve. 

The total diode current density is then 

00 „ 

I tot = (eD/irh) £ / T L ^ R (E 0x ,L ) ■ f FD (E 0x , L , fi(B))dE 0x -[£<-> R], (16) 

L =0 J 
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where the second term in the right-hand side is the backward contribution to the total 
current. In understanding the abnormal characteristics of the current, it is helpful to define 
an effective chemical potential fi e ff(B, Lq) dependent upon the Landau level Lq as follows 



At extremely low temperatures such as 4.2 K, the Fermi-Dirac function becomes flat for 
energies below the chemical potential. To contribute to the total current, the longitudinal 
energy of electron states (E 0x ,L ) should satisfy E 0x < fi e ff(B,L ). 

As is indicated by Eq.(16) the resonant tunneling current through a double barrier struc- 
ture generally depends on the details of the transmission probability, the density of electronic 
states (DOS) and the chemical potential fi (or Fermi energy). Among these, the DOS is much 
more affected by strong magnetic fields rather than the other properties. As the summation 
and factor 1 / y/ E 0x in Eq.(14) indicate, the DOS severely changes due to the quantization 
of the lateral energy of electrons in the left (right) local contact. The electron states are 
then classified using the longitudinal sub-band index Lq. At the temperature of 4.2 K, the 
longitudinal energy range depends on L Q according to E 0x < fi e ff(B, Lq). Thus the effective 
chemical potential fi e ff(B, L ) for a particular sub-band L is defined to indicate the cutoff 
obtained in the longitudinal energy space. 

Under such conditions, the current for an occupied sub-band L reduces to 



where the minimum of the conduction band at the left contact has been set to zero. There- 
fore, the area of the transmission probability determines the contribution of each occupied 
sub-band to the total current. Fig. shows the currents calculated for structure A and B. 
As the voltage is increased, the lowest L sub-band contributes to the current first with 
higher Landau levels contributing at higher voltages. 



fj, eff (B, L ) = n{B) - (2Lq + l)h Wd0 /2 



(17) 




T L ^ R (E 0x , L )dE 0x - [L^R] 



(18) 
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III. RESULTS 



We note here two important features in the results obtained. First, a shift of the onset 
voltage depending on the magnetic field strength is observedS. It is sharper for the structure 
with thicker barriers. Second, the current peaks are enhanced in regions of strong fields in 
both structures^. For a thicker barrier (structure B) this current increase is quite remarkable 
over a wide range of B field. 

Let us analyze why these features occur in the presence of a magnetic field. As the applied 
voltage increases, the quasi-resonant energy E r (see the inset of Fig. p is lowered. At the 
same time, the conduction band profile has a decreased symmetry due to its tilt. It results 
that the height and width of the peak transmission probability near E r decreases, as can 
be verified in Fig. [I]. The energy width of the transmission probability is narrower in RTD 
structures with a thicker barrier. Therefore the transmission probability is much sharper for 
structure B than for structure A. These properties explain the increased sharpness of the 
current threshold and the downward decline of the current plateau in structure B (see Fig. 

Let us consider the lowest Landau sub-band Lq = 0, whose effective chemical potential 
is fi e ff(B, 0) = /i(-B) — hw c0 (B)/2. In the absence of a magnetic field, the onset voltage for 
RTD diodes with symmetric barriers and for eVp = Ef^ —Ef ^ positive is well-known to be 
approximately given by 

V D)0n /2 = E r - E f)L = E r - MO), (19) 

where Efi is the Fermi energy level for the electron carriers on the left contact. In the 
presence of a magnetic field /z(0) is just replaced by /i e //(i?,0) and E r can be assumed to 
remain unchanged. Indeed the dependence of E r upon the magnetic field B which arises 
from the variation of the mass along the superlattice direction (see Eq.(8)) was verified to 
be negligible in the structures studied. As a result the shift AVd i0 „ of the onset voltage is 
approximately given by 
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AVb )0n /2 = [V D>on {B) - Vd, o „(0)]/2 ~ [/i(0) - fx(B)] + hw c0 (B)/2, (20) 

which is equal to the difference in effective chemical potential. Fig. ^ indicates that the 
dominant contribution to the onset voltage shift originates from the non-zero energy hwco/2 
of the Landau ground (lowest) state which is a unique feature. Note also that the superposed 
oscillation in AVr> y0n has the same origin, namely the oscillation of the chemical potential, 
as in other phenomena such as de Haas-van Alphen 

Finally let us look at the increase of the current peak (see Fig. |2]). The current density 
Eq.(18) has a weight factor proportional to the degeneracy D. As discussed, the sharpness of 
the transmission coefficient depends on the structure parameters such as the barrier widths 
as well as the conduction-band barrier heights. The enhancement of the current peaks in 
presence of magnetic field depends on the width of the transmission coefficient peak near the 
quasi- resonant energy E r . Therefore, a structure with a sharper transmission probability 
has a more enhanced peak. 

The RTD diode current in Fig. |2| for structure B increases step by step for a weaker 
magnetic field, with plateaus slowly declining downward. However for a strong magnetic 
field at B = 20 Telsa where the two lowest Landau levels of the emitter are occupied, its 
diode current rises to its maximum peak value in just one step. For structure A (B) in a 
strong field of B = 20 Telsa, an enhancement of the current peak by 10(37) % is observed. 



IV. CONCLUSION 

In this paper we have analyzed coherent electron transport in double-barrier heterostruc- 
tures with parallel electric and magnetic fields with the aid of a quantum simulator 
accounting for 3-dimensional transport. The presence of a magnetic field affect severely 
the electron transport in a resonant tunneling diode causing a shift of the onset voltage. 

Our analysis establishes that the onset-voltage shift induced by the magnetic field arises 
from an upward shift of the non-zero ground (lowest) Landau state energy in the entire 
quantum region where coherent transport takes place. Note that quantum region therefore 
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includes the left contact rather than only the well as previously reportedu. Actually in this 
work the left contact is used as the cyclotron energy reference. Indeed our 3-dimensional 
transport analysis accounts for the spatial dependence of the cyclotron frequency which 
is, however, verified to have a relatively small impact on resonant tunneling in the device 
analyzed for the magnetic field strength considered. A larger correction term for the onset- 
voltage shift arising from the magnetic field dependence of the chemical potential is also 
derived. 

The comprehensive 3D analysis reported in this paper reveals that the Landau ground 
state with its nonvanishing finite harmonic oscillator energy hw c (0)/2 is indeed the principal 
contributor to the onset voltage shift at low temperatures. Finally it is also shown in this 
paper how the abnormal increase of the current peak depends on the external magnetic field 
as well as the structure parameters of the double barrier heterostructure. 
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FIGURES 

FIG. 1. The inset shows the conduction band profile of the resonant-tunneling structure. The 
transmission probability with increased voltage at T = 4.2K, magnetic field B = 10 Telsa, and the 
Landau state Lq = for structure A and B. The solid, dotted, and dashed lines are calculated for 
V D = 0.18,0.23,0.28 V for structure A and V D = 0.07,0.12,0.17 V for structure B, respectively. 

FIG. 2. The current-voltage characteristic at T = 4.2K for the structure A and B in pres- 
ence of magnetic fields. The solid, dotted, dashed, dot-dashed, double dot-dashed lines are for 
B = 0,5, 10, 15, 20 Telsa, respectively, for both structures. 

FIG. 3. The solid line is the shift of onset voltage at T = 4.2K in structure B for increasing 
magnetic field. The dashed line shows the contribution in Eq.(20) of the energy l/2hw c of the 
Landau ground (lowest) state. 
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